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ABSTRACT 



We present a spectroscopic search for halo field stars that originally formed in globular clusters. Using moderate-resolution SDSS- 
III/SEGUE-2 spectra of 561 red giants with typical halo metallicities (-1.8 <[Fe/H] < -1.0), we identify 16 stars, 3% of the sample, 
with CN and CH bandstrength behavior indicating depleted carbon and enhanced nitrogen abundances relative to the rest of the data 
set. Since globular clusters are the only environment known in which stars form with this pattern of atypical light-element abundances, 
we claim that these stars are second-generation globular cluster stars that have been lost to the halo field via normal cluster mass- 
loss processes. Extrapolating from theoretical models of two-generation globular cluster formation, this result suggests that globular 
clusters contributed significant numbers of stars to the construction of the Galactic halo: we calculate that a minimum of 17% of the 
present-day mass of the stellar halo was originally formed in globular clusters. The ratio of CN-strong to CN-normal stars drops with 
Galactocentric distance, suggesting that the inner-halo population may be the primary repository of these stars. 

Key words. Stars: abundances - Galaxy: halo - Galaxy: formation 



1. Introduction 

Light-element abundance inhomogeneities in globular clusters 
have been a subject of ongoing study for over thirty years. Early 
observations of a bimodal distribution of CN bandstrength in 
globular cluster red-giant (RGB) stars, antic orrelated with CH 
bandstrength (e.g., Norris & Cottrell 1979; Suntzeff 198 1 ; Norris 
et al. 1984), were interpreted as evidence of anticorrelated ranges 
of carbon and nitrogen abundances in globular cluster stars (e.g., 
Bell & Dickens 1980). Since that time, observations have been 
extended to include more clusters, fainter stars within individ- 
ual clusters, and additional elemental abundances. We now know 
that anticorrelated C-N, O-Na and Mg-Al ranges are present in 
stars on the red-giant branches of all globular clusters that have 
been thoroughly studied (e.g., Kayser et al. 2008; Carretta et al. 
2009; Smolinski et al. 201 lb), and that the C-N anticorrelation is 
also found among main-sequence stars in several clusters (e.g., 
Harbeck et al. 2003; Briley et al. 2004; Pancino et al. 2010; 
Smolinski et al. 2011b). Light-element abundances in globular 
cluster stars range from scaled-Solar, similar to halo field stars of 
the same metallicity, to depleted in C, O and Mg and enha nced 
in N, Na and Al, with be tween 30% dPancino et al.ll2010l) and 
70% (ICarretta et al.ll2009l) of stars in a given cluster having an 
atypical abundance pattern. 

The atypical abundance pattern resembles the result of high- 
temperature hydrogen burning: the CNO cycle converts carbon 



and oxygen into nitrogen, and the NeNa and MgAl cycles in- 
crease the abundances of sodium and aluminium while deplet- 
ing magnesium. The presence of this abundance pattern in old, 
low-mass main-sequence stars, which have neither the high- 
temperature fusion zones nor the capacity to transport material 
between their cores and surfaces, requires that the abundance 
variations be extrinsic. The fact that the abundance differences 
between stars with typical and atypical abundance patterns are 
equally large before and after first dredge-up, which occurs low 
on the RGB, indicates that the abundance differences are present 
throughout the stars, and are not merely surface pollution. 

As a result, the prevailing model for the origin of primor- 
dial light-element abundance variations in globular clusters is 
that the stars with atypical light-element abundances are a nearly 
coeval second generation formed from material processed by 
intermediate- or high-mass stars in the first generation. There 
are several possible sites for high-temperature fusion process- 
ing in the first generation. AGB stars are a common suggestion 
(e.g., Parmentier et al. 1999) because they have relatively slow 
winds, they are a site of hot hydrogen burning (e.g., Karakas 
2010), and they evolve on timescales of ~ 10 s years, quite fast 
compared to the lifetime of a globular cluster. Other possible 
sources for high-temperature- processed feedback m aterial are 
rapidly rotating massive stars dDecressin et al.ll2007l) and mas - 
sive binary stars undergoing mass transfer (Ide Mink et alj2 009). 
both of which could return more feedback mass in a shorter 
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Fig. 1. Surface gravity versus effective temperature for our initial 
(small gray points) and trimmed (open gray triangles) data sets. 
Solid black lines are 12 Gyr Dartmouth isochrones with (left to 
right) [Fe/H]= -1.8, -1.4, and -1.0. 



amount of time than AGB stars, although the stars themselves 
are less common. 

However, all of these possible feedback sources suffer from 
what has come to be known as the "mass budget problem," 
the mismatch between the current 1 : 1 ratio of first- to second- 
generation stars and the concept that the second generation was 
built entirely from winds from the first generation. Recent the- 
oretical models for globular cluster formation propose two so- 
lutions to the mass budget problem: massi ve first generations 
and s ignificant gas a ccret ion. The models of Be kki et all (2007) 
and iD'Ercole et al.1 (120081) both require a first generation with 
an initial mass 10 to 20 times its present mass, and have sec- 
on d generations built en tirely from AGB winds. In the model 
of ID'Ercole et al.l (|2008), the second-generation star formation 
occurs near the center of the cluster, and concludes when type 
la supernovae begin. The la SNe drive an overall expansion of 
the cluster, causing stars at large radii (primarily first-generation 
stars) to dissociate from the cluster, thereby reducing the ratio of 
first- to second-generation stars to present-day levels. AGB stars 
are also the source of chemical differen ces between the firs t and 
second generation in the model of lConrov & Sperg"efl(l201 lb . but 
their model involves a significant amount of gas being accreted 
from the environment, reducing the mass required for the first 
generation of stars. The authors investigate the mass-accretion 
rate for proto-globular clusters in a variety of cosmological en- 
vironments, and conclude that clusters with initial masses above 
1O 4 M would have been able to accrete significant amounts of 
gas in the early Milky Way. They also calculate that the mini- 
mum mass for sufficient gas accretion would be lower for glob- 
ular clusters forming as satellites of dwarf galaxies, a claim that 
could be tested with observations of intermediate-age star clus- 
ters in the Large Magellanic Cloud (LMC). Inte r esting ly, the re- 
cent high-resolution study of iMucciarelli et ail (1201 ll) finds no 
sign of light-element anticorrelations among 14 red giant stars 
in the massive intermediate-age LMC cluster NGC 1866. The 
authors suggest that the minimum mass for self-enrichment in 
the LMC is on the order of 10 5 M o . Although there are ques- 



Fig. 2. Two sample spectra from our data set, demonstrating the 
effects of varying [C/Fe] and [N/Fe] at fixed stellar parame- 
ters. The upper spectrum, of SDSS J1500+1 126, is CN-strong 
while the lower spectrum, of SDSS J1524+0550, is CN-normal. 
Shaded regions denote the bandpasses of 5(3839) and S{CH): 
lines angled up to the right mark the sideband (more closely 
spaced lines) and science band (more broadly spaced lines) of 
5(3839), and lines angled down to the right mark the sidebands 
(more closely spaced lines) and science band (more broadly 
spaced lines) of S(CH). 



tions about the likelihood that early globular clusters would ac- 
crete unprocessed gas from their surroundings with precisely the 
right metallicity (e.g., Martell 2011), and the gas dynamics in- 
volved in mass loss and dilution in globular clusters have been 
the subject of few advanced numeri cal studies (e.g., Priestley 
et al. 2011; D'Ercole et al. 2011), the lConrov & Spergell (1201 11) 
model provides an important discussion of the interplay between 
early cluster self-enrichment and the galaxy-scale environment. 

These scenarios rely implicitly on the high mass and high 
density of early globular clusters to retain and/or accrete gas well 
enough to permit a second burst of star formation. It is there- 
fore not surprising that although stars with second-generation- 
like light-element abundances are found in every globular clus- 
ter in the Milky Way, they are not found at all in old open clus- 
ters. In studies of individual stars (e.g., Jacobson et al. 2008; 
Martell & Smith 2009) and mean abundances (e.g., de Silva et al. 
2009) for old open clusters, light-element abundance behavior is 
found to be distinctly different from what is observed in globular 
clusters. They are also quite unlikely to have formed in the halo 
field, b ut they have recently been found there: Martel T& Grebell 
d2010h (hereafter MG10) searched the SEGUE survey (Sloan 
Extension for Galactic Understanding and Exploration, Yanny 
et al. 2009) for halo giants with unusually strong UV/blue CN 
bands and identified 49 (of roughly 2000) stars likely to have 
low carbon abundanc es and high nitrogen abundances. Similarly, 
ICarretta et al.l d2010l) compiled a sample of 144 metal -poor disk, 
halo and bulge stars from the literature and identified 2 of those 
as Na-rich and likely to have originated in globular clusters. 
In both cases, the authors interpret the field stars with second- 
generation abundances as stars that formed in globular clusters 
and were later transferred to the halo through mass-loss pro- 
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Fig. 3. CN bandstrength index 5(3839) versus absolute r mag- 
nitude for our full data set, divided into 0.1-dex-wide bins in 
[Fe/H]. Maximum [Fe/H] in each bin is given in the upper left 
corner of each panel. Solid lines are best fits to the CN-normal 
stars in each panel. The slope of the best-fit line and the separa- 
tion between CN-normal and CN-strong stars both decline with 
decreasing metallicity. 



cesses such as tidal stripping by the Galaxy or two-body in- 
teractions within the cluster. It is also claimed in MG10 that 
their figure of 2 .5 % of halo stars having second-generation abun- 
dances, seen through the lens of a globular cluster formation 
model with strong early mass loss (such as, e.g., D'Ercole et 
al. 2008), implies that as much as 50% of halo field s tars origi- 
nally formed in globular clusters. ICarretta et al.l (T2010) conclude 
that early mass loss from globular clusters is "a major building 
block of the halo", and may also have contributed significantly 
to the formation of the thick disk. This possibility has been sug- 
gested theoretica l ly as well: both iBaumgardt et al.l d2008l) and 
Mark s & Kroupal Q2010) find that early mass loss from globular 
clusters should contribute significant numbers of stars to the halo 
field. 

In this paper we use data from the SEGUE-2 moderate- 
resolution spectroscopic survey, available as part of the eighth 
data release of the Sloan Digital Sky Survey (Aihara et al. 201 1; 
Eisenstein et al. 201 1; York et al. 2000), to search for CN-strong 
halo giants like those identified in MG10. This will allow us to 
make two important extensions to that work: first, by expand- 
ing the sample of stars surveyed, we can refine our estimate of 
f^ G , the present-day fraction of CN-strong stars in the halo field; 
second, we calculate f GC , the fraction of halo stars with first- or 
second-generation chemistry originating in globular clusters, as 
predicted by several different two-generation cluster formation 
models. 



2. The Data Set 

The SDSS and its extensions have acquired ugriz pho- 
tometry for several hundred million stars. SEGUE, one of 
three sub-surveys that together formed SD SS-II, extended 
the ugriz imaging footprint of S DSS-I dFukugita et alj 1 1996; 
iGunn et all [1998H York et alll2000l: IPier et al.ll2003b iGunn et alJ 



Fig. 4. CN-CH planes for our final data set, divided into 0. 1-dex- 
wide bins in [Fe/H]. The CH bandstrength index S(CH) is plot- 
ted versus the CN bandstrength residual 6S (3839), and candidate 
globular cluster stars are drawn as open circles. As in Figure 3, 
the maximum metallicity in each bin is given in the upper left 
corner. As is explained in the text, we select as candidate glob- 
ular cluster stars all those with relatively high 55(3839) in the 
three highest-metallicity bins, and stars with both high (55(3839) 
and S (CH) below the mean S (CH) of the CN-normal group in 
all other metallicity bins. 



2006; Stought on et al.ll2002t lAbazaiian et al.ll2003ll2004 120051 



Adelman-McCarth v et alJ I2006L l2007l 120081: lAbazaiian et al.l 
2009) by approximately 3500 deg 2 , and also obtained R 2000 
spectroscopy for approximately 240, 000 stars over a wavelength 
range of 3800 - 9200 A. SEGUE included spectra for a collection 
of Galactic globular and open clusters, which served as calibra- 
tors for the T e ff, log(g), and [Fe/H] scales for all stars observed 
by SDSS/SEGUE, as processed by the SEGUE Stellar Parameter 
Pipeline (SSPP; Lee et al. 2008a; b; Allende Prieto et al. 2008). 
The SSPP produces estimates of T e ff, log(g), [Fe/H], and radial 
velocity, along with the equivalent widths and/or line indices 
for 85 atomic and molecular absorption lines, by processing 
the calibrated spectra gen erated by the standard SDS S spectro- 
scopic r eduction pipeline (Stough ton et al.l l2002). See lLee et al.l 
(2008a) f or a detailed discuss ion of the approaches used by the 
SSPP: ISmolinski et al.l (1201 1 ah provides details on the most re- 
cent updates to this pipeline, along with additional validations. 
SEGUE-2 (C. Rockosi et al., in prep.) expanded the numbers 
of Galactic stars with available low-resolution SDSS spectra by 
over 120,000, concentrating in particular on stars with greater 
distances in order to better sample the outer-halo region of the 
Galaxy. 

As in MG10, SEGUE-2 targets were selected from across a 
large range in parameter space: we began by selecting all SDSS- 
III targets from SEGUE-2 plates with [Fe/H]< -1.0, log(g) < 
3.0, (g - r)o > 0.2, and a mean signal-to-noise ratio (SNR) of at 
least 20 per pixel. Errors on the derived quantities and quality- 
assurance flags were also required to fall within a certain range: 
o"[Fe/H] ^ 0.5, CH g( g ) < 0.5, cr Teff < 200 K. This initial data set 
contained 2019 stars. We then identified four stars as carbon-rich 
using the C2 indices defined in MG10, and removed them from 
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the sample. RGB stars were identified by dividing the data into 
0.2-dex-wide metallicity bins and making a recursive 3-sigma 
selection about a fiducial sequence in the (log(g), T e ff) plane. 
Since our analysis is based on the CN and CH molecular ab- 
sorption features, which become quite weak at low metallicity, 
we also restrict our sample to stars with [Fe/H] > —1.8 and a 
SNR per pixel between 4000 A and 4 100 A (SN Wue ) above 15. 
The RGB, metallicity, and SNbi ue selections together reduce our 
data set to 561 stars. 

Some stars from SEGUE- 1 were intentionally re-observed 
in SEGUE-2, as a way to check for consistency in data reduc- 
tion and analysis between the two surveys. However, none of the 
stars in this study were also in the MG10 data set. Table 1 lists 
SDSS identifiers, the plug plate number, MJD, and fiber number 
of the SEGUE observation (useful for selecting individual stars 
from the SDSS Catalog Archive ServetQ), right ascension, decli- 
nation, photometry, SSPP-derived stellar parameters, and survey 
name (SEGUE- 1 or SEGUE-2) for each star in the final trimmed 
data set and the final data set from MG10. Only the first ten 
rows of Table 1 are given in the print version of this paper, to 
indicate form and content; a full version is available at the CDS 
0. Figure 1 shows surface gravity versus effective temperature 
for our initial (small gray points) and trimmed (open gray tri- 
angles) data sets, using the T e fr and log(g) values derived from 
the spectra by the SSPR The overplotted solid lines are 12 Gyr 
Dartmouth isochrones (Pott er et alJl2008h with metallicities of 
[Fe/H]= -1.8, -1.4, and -1.0, spanning the metallicity range of 
our data set. 

3. Analysis 

Because we are searching for old, low-mass red giant stars with 
globular-cluster-like abundance patterns, our analysis follows a 
process typical in the globular cluster literature. Specifically, we 
are using index-based techniques developed for globular clus- 
ter studies as indicators of carbon and nitrogen abundance, and 
we infer the presence of the full second-generation light-element 
abundance pattern from those indices. Overall, we followed the 
procedure from MG10 fairly closely in analyzing our data, in 
order for the results to be as compatible as possible. 

3.1. Bandstrength indices 

We use t he bandstrength ind ices 5 (3839) (iNorris et alJl98ll) and 
S(CH) (Marte ll et al.ll2008l) . which were initially developed for 
studies of light-element abundance variations in globular clus- 
ters, to identify stars with strong CN bands and weak CH bands 
relative to the majority of the halo field. Indices measure the 
magnitude difference between the integrated flux in the feature 
(the "science band") and one or two nearby, independent regions 
of the spectrum (the "sidebands"), in the sense that more absorp- 
tion in the feature results in a larger index value. Since molec- 
ular abundance is strongly controlled by the abundance of the 
minority species, the CH band is indicative of [C/Fe] while the 
CN band traces [N/Fe]. We are using carbon and nitrogen abun- 
dances as an indicator of the "second-generation" abundance 
pattern of depleted carbon, oxygen, and magnesium and en- 
hanced nitrogen, sodium, and aluminum because that full abun- 
dance pattern is consistently found together in Galactic globular 
clusters. 

1 http://skyservice.pha.jhu.edu/CasJobs/login.aspx 

2 via anonymous ftp to cdsarc.u-strasbg.fr (130.79.128.5) or via 
http://cdsweb.u-strasbg.fr/cgi-bin/qcat?J/A+A/ 



Figure 2 shows typical spectra from our data set, chosen to 
illustrate the range of variation in CN and CH bandstrength. 
The upper spectrum (of SDSS J1500+1 126) is of a CN-strong 
star with T eff = 5250, [Fe/H]= -1.70, 5(3839) = 0.005 and 
S(CH) = 0.70, and the lower spectrum (of SDSS J1524+0550) 
has similar stellar parameters but relatively weaker CN and 
stronger CH bands (T eff = 5210, [Fe/H]= -1.74, 5(3839) = 
-0.189 and S(CH) = 0.79). The spectra are both normalized 
near 4500A and offset vertically. The shaded areas with lines 
tipping up to the right mark the sideband (more closely spaced 
lines) and science band (more broadly spaced lines) of 5(3839), 
and the shaded areas with lines angled down to the right mark the 
sidebands (more closely spaced lines) and science band (more 
broadly spaced lines) of S(CH). The bandpasses of both indices 
are given in Table 2. Errors on the bandstrengths were estimated 
as in MG10, using a Monte Carlo technique to randomly sample 
the SEGUE noise vectors and add that noise to the spectra, then 
measure the standard deviation of the bandstrength over 100 re- 
alizations. Typical values for os^g) and crs(CH) are 0.015 and 
0.007, respectively. 



3.2. Identification of CN-Strong candidates 

We identify candidate second-generation halo stars in our data 
set in two stages, first by selecting those stars with unusu- 
ally high 5(3839) relative to other stars of similar metallicity 
and evolutionary stage, and then by taking just the subset of 
those CN-strong stars with appropriate S{CH). Absolute mag- 
nitudes for our stars were calculated using a spectroscopic par- 
allax method similar t o that used in MG 10. We used 12-Gyr 
Dartmouth isochrones dDotter et al.l 1200 8), interpolated to the 
metallicities of our target stars, then use dereddened (g - r)o col- 
ors to predict absolute M ( . magnitudes. 

Figure 3 shows 5 (3839) versus M r for our 56 1 program stars, 
divided into 0.1-dex-wide metallicity bins. The solid lines over- 
plotted are linear fits to the CN-normal stars, and the maximum 
metallicity in each bin is given in the upper right corner of each 
panel. The distinctive characteristics of each panel are the dom- 
inance of the CN-normal group, the slope of the best-fit line, 
and the separation of the CN-normal and CN-strong groups. In 
globular clus ters, the fraction of CN-strong stars is roughly 50% 
(Ki-afi ll994h . but the slope of the CN bandstrength-luminosity 
relation and the separation in CN bandstrength between the two 
groups behaves similarly to what we observe in field stars. As in 
MG10, the slope and the separation both decline with decreasing 
metallicity, and the two CN bandstrength groups become diffi- 
cult to distinguish at the lowest metallicites. 

We measure the quantity (55(3839), the difference between 
an individual star's CN bandstrength and the best-fit line at the 
same absolute magnitude, in the plane of Figure 3. The eight 
panels of Figure 4 correspond to the panels of Figure 3, and 
show <55(3839) versus S(CH) for each metallicity bin, with can- 
didate second-generation field stars plotted as open circles and 
all other stars plotted as filled circles. Because the S(CH) in- 
dex has a reduced sensitivity to carbon abundance at the upper 
end of our metallicity range, we accept all CN-strong stars with 
metallicities above [Fe/H] = -1.3 as candidates, and we select 
only those CN-strong stars with 5 (CH) below the mean of the 
CN-normal group for stars with metallicity below [Fe/H] = -1.3. 
These are the same criteria as were used in MG10, and they re- 
turn 16 candidate globular cluster stars, 3% of our sample. This 
is very similar to the MG10 result based on SEGUE-I stars. 
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Table 1. Star identifiers, position, photometry and parameters for all program stars from this study and from MG10. 
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SDSS ID 


Plate 


MJD 


FiberlD 


a 


6 


£o 


(g ~ r) 


[Fe/H] 




log(g) 


Survey 


SDSS J1605+1638 


3289 


54910 


132 


16:05:52.72 


16:38:59.21 


16.04 


0.663 


-1.07 


4965 


2.11 


SEGUE-2 


SDSS J0842+2216 


3373 


54940 


51 


08:42:16.23 


22:16:58.27 


16.14 


0.715 


-1.05 


4913 


2.06 


SEGUE-2 


SDSS J1528+0759 


3308 


54919 


529 


15:28:51.34 


07:59:08.62 


17.11 


0.696 


-1.04 


4950 


2.29 


SEGUE-2 


SDSS J1557+3629 


3478 


55008 


229 


15:57:48.95 


36:29:33.04 


16.77 


0.735 


-1.03 


4910 


2.29 


SEGUE-2 


SDSS J1011+1810 


3178 


54848 


269 


10:11:13.50 


18:10:32.73 


17.22 


0.689 


-1.00 


4934 


2.37 


SEGUE-2 


SDSS J2243+1322 


3128 


54776 


104 


22:43:31.49 


13:22:33.68 


17.72 


0.689 


-1.09 


4877 


2.55 


SEGUE-2 


SDSS J1016+1737 


3178 


54848 


55 


10:16:42.76 


17:37:23.51 


17.44 


0.619 


-1.05 


5066 


2.50 


SEGUE-2 


SDSS J1206+3021 


3181 


54860 


73 


12:06:49.60 


30:21:49.39 


18.06 


0.739 


-1.05 


4983 


2.53 


SEGUE-2 


SDSS J1210+1042 


3214 


54866 


14 


12:10:11.17 


10:42:28.64 


16.80 


0.652 


-1.07 


5025 


2.56 


SEGUE-2 


SDSS J1617+0502 


3298 


54924 


48 


16:17:42.83 


05:02:32.70 


17.54 


0.663 


-1.06 


4991 


2.56 


SEGUE-2 
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Fig. 5. Heliocentric distance as a function of absolute r magni- 
tude for our final data set (upper panel) and for the MG 1 data set 
(lower panel). CN-normal stars are plotted as small dots, while 
CN-strong stars are shown as filled stars. Only a randomly se- 
lected 25% of the CN-normal MG10 stars is shown to reduce 
crowding. In both panels, the lower boundary is set by the bright 
limit of the SEGUE survey (g = 14), and the upper boundary is 
set by our requirement that SNbiue > 15. 



3.3. Stellar distances 

We calculate heliocentric distances D to our target stars from the 
distance modulus (rrj— M r ). Errors in distance were calculated by 
randomly sampling the error on (g - r)o, then recalculating the 
distance. The standard deviation in 100 of these realizations was 
then adopted as <x D . One of the goals of the SEGUE-2 survey 
was to observe stars at greater distances than were observed in 
SEGUE-1, and while our data set covers a heliocentric distance 
range very similar to that of MG10, the average star at any fixed 
M r is more distant for the current data set than for the MG10 data 
set. The upper panel of Figure 5 shows our calculated distances 
versus absolute M r magnitudes for the SEGUE-2 data set, with 
CN-normal stars drawn as small dots and CN-strong stars drawn 
as filled stars. The lower panel shows the analogous data for the 
MG10 data set, with only a randomly selected 25% of the MG10 
CN-normal stars plotted for visual clarity. In both panels, the 
stars fill a band that is restricted at small distances by the SEGUE 



Fig. 6. CN bandstrength index 55(3839) versus Galactocentric 
distance roc for our final data set (upper panel) and the MG10 
sample (lower panel). Symbols are the same as in Figure 5. There 
is a sharp drop in the number and frequency of CN-strong stars 
at roughly 20 kpc. 



bright limit of g — 14 and at large distances by our requirement 
that SN b i U e > 15.0. 

Heliocentric distances are converted geometrically 
to Galactocentric distances roc, using the relation 
r GC = D 2 + (8kpc) 2 - 2D cos(y), in which y is the angle 
on the sky between the Galactic center and the star in question. 
Table 3 lists plate, MJD, and FiberlD identifiers, along with 
our measurements of S (3839) and S (CH), our determination of 
a star's CN class (normal or strong), absolute M r magnitude, 
heliocentric distance, Galactocentric distance, and survey name 
(SEGUE-1 or SEGUE-2), for the same stars as in Table 1. As 
with Table 1 , only the first ten rows are given in the print version 
of this paper, and a full version of the table is available from 
CDS. Figure 6 shows 55(3839) versus rcc for our final data set 
(upper panel) and the MG10 sample (lower panel), using the 
same symbols as in Figure 5. The frequency of CN-strong stars 
appears to drop at roughly 20 kpc, with only one CN-strong star 
beyond rcc = 30 kpc. 

To further investigate this result, we visualize the re- 
lationship between the frequency of CN-strong stars and 
Galactocentric distance in two ways: Figure 7 shows f s , the ra- 
tio of CN-strong to CN-normal stars, versus Galactocentric dis- 
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Fig. 7. The ratio of CN-strong to CN-normal stars, f„ as a func- 
tion of Galactocentric distance for the present data set and the 
MG10 data combined. Vertical error bars represent the Poisson 
errors on f s in each distance bin. The drop in the frequency of 
CN-strong stars seen in Figure 5 at 20 kpc is also seen here. 



Fig. 8. Cumulative distribution functions for (55(3839) for stars 
with rcc < 20 kpc (solid line) and rcc ^ 20 kpc. The fraction of 
stars rises more quickly at large values of 6S (3839) for stars with 
smaller Galactocentric radii, indicating that CN-strong stars are 
more common in the inner halo. 



Table 2. Bandstrength index definitions 



Index 


Blue sideband 


Science band 


Red sideband 


5(3839) 
S(CH) 


4050-4 100 A 


3846-3883A 
4280-4320A 


3883-3916A 
4330-4350A 



tance for the present data set combined with MG10. Vertical er- 
ror bars represent Poisson errors in f s . The fraction of CN-strong 
stars is lower at larger rcc, with a transition at a 20 kpc: rela- 
tively small-rcc stars are consistent with a value for f s around 
0.04 and relatively large-rcc stars are consistent with a value 
for f s around 0.01. Figure 8 shows cumulative distribution func- 
tions of 55(3839) for stars with rcc < 20 kpc (solid line) and 
rcc ^ 20 kpc (dotted line). The cumulative distribution curve 
for stars with smaller Galactocentric distances climbs continu- 
ously for 55(3839) > 0.05, while the curve for more distant 
stars is flatter in that range, indicating that there are relatively 
more CN-strong stars at smaller Galactocentric distances. A one- 
sided Kolmogorov-Smirnovtest returns a probability of 0.15 that 
the 55(3839) distributions in these two distance ranges were 
drawn from the same parent population. This result is only a very 
marginal rejection of the single population hypothesis; however, 
it clearly points to a reasonable possibility that two parent popu- 
lations with differing CN-strong fractions exist. 

We interpret this drop-off in the frequency of CN-strong 
stars with Galactocentric distance as a possible sign of a tran- 
sition from the inner-halo population to the outer-halo popu- 
lation (as discussed by, e.g., Carollo et al. 2007; 2010). Note 
that this also corre sponds to the transition zone reported by 
Ide Jong et alJ d2010h . based on fits to Hess diagrams from the 
SEGUE vertical photometry stripes. Recent theoretical studies 
of galaxy formation (e.g., Oser et al. 2010; Font et al. 201 1) pre- 
dict that stars formed in situ are the dominant population at small 
Galactocentric radii, while stars accreted during minor mergers 
are the majority at larger Galactocentric distances. Studies of 
abundances in individual stars in dwarf galaxies do not often 



report stars with light-element abu ndances similar to se cond- 
generation globular cluster stars. In lShetrone et al.l d2003) there 
is only one star out of the 15 surveyed (Forn ax 21) with even 
mildly elevated [Na/Fe], Sadakan e et al.l d2004t) found fairly low 
[Na/Fe] ab undances in t h e thre e stars they observed in Ursa 
Minor, and iGeisler et al. (120051) find the same result for nine 
RGB stars in Sculptor. Koch et al] d2008l) also report one star 
in Carina (of eight wi th measured sodium abundance) with 
moderate [Na/Fe], and iLetarte et al.l (1201 Oh report one star in 
Fornax (of 70 with measured sodium abundances) with mod- 
erate [Na/Fe], but neither study finds any stars with the level 
of sodium abundance enhancement typically found in second- 
generation globular cluster stars (e.g., Carretta et al. 2009). 

A schematic picture of the Galaxy in which the inner halo 
has significant contributions from globular clusters, while the 
outer halo is dominated by stars that formed in dwarf galaxies, 
fits reasonably well with our observation that a higher fraction of 
stars in the inner halo have light-element abundances similar to 
those of second-generation globular cluster stars. Spectroscopic 
surveys currently underway (such as APOGEE; Allende Prieto 
et al. 2008), surveys presently being planned as companions to 
the Gaia mission (e.g., Cacciari et al. 2009), and studies based 
on nearby stars with accurate proper motions and halo kinemat- 
ics (such as, e.g., Carollo et al. 2010) will be important ways to 
explore this question of inner/outer halo origins further: all of 
them expand the effective search space, either by surveying to 
larger distances or by focusing on stars with large apogalactic 
distances. 



4. Interpretation 

The homogeneity of metallicity in present-day globular clusters 
seems to favor massive models of globular cluster formation, 
while the observed range in oxygen and sodium abundances 
(see, e.g., Carretta et al. 2009) indicates that dilution of feed- 
back material and incomplete mixing of the gas that makes up 
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Table 3. Measured and derived parameters for all program stars from this study and from MG10. 



7 



Plate MJD FiberlD 5(3839) S(CH) CN class M r D (kpc) R GC (kpc) Survey 



3289 


54910 


132 


-0.026 


0.716 


Normal 


0.31 


10.3 


17.7 


SEGUE-2 


3373 


54940 


51 


0.004 


0.760 


Normal 


-0.12 


12.9 


17.2 


SEGUE-2 


3308 


54919 


529 


-0.008 


0.783 


Normal 


0.06 


18.6 


25.2 


SEGUE-2 


3478 


55008 


229 


0.014 


0.754 


Normal 


-0.23 


18.0 


25.5 


SEGUE-2 


3178 


54848 


269 


0.035 


0.763 


Normal 


0.22 


18.2 


22.6 


SEGUE-2 


3128 


54776 


104 


0.025 


0.710 


Normal 


0.01 


25.3 


30.3 


SEGUE-2 


3178 


54848 


55 


-0.153 


0.789 


Normal 


0.80 


15.9 


20.5 


SEGUE-2 


3181 


54860 


73 


-0.106 


0.819 


Normal 


-0.32 


33.8 


38.2 


SEGUE-2 


3214 


54866 


14 


-0.083 


0.798 


Normal 


0.40 


14.1 


19.3 


SEGUE-2 


3298 


54924 


48 


-0.105 


0.784 


Normal 


0.33 


20.4 


27.6 


SEGUE-2 



the second stellar generation play a central role in setting final 
light-element abundances for individual stars. In order to evalu- 
ate the veracity of various globular cluster formation models, it 
is necessary to explore their implications for the formation of the 
Galactic stellar halo. Massive models, which require a large frac- 
tion of first-generation stars to be lost to the halo field at early 
times, will produce a lar ge enhancement between fl G a nd f GC . 

The recent paper of Scha erer & Charbonnelld201ll) consid- 
ers this question from a theoretical sta ndpoint, using th e glob- 
ular cluster self-enrichment model of iDecressin et al.1 (|2007) 
to calculate the mass in first-generation stars needed to pro- 
duce a significant second generation. In the IDecressin et al.l 
(2007) model, rapidly rotating massive stars provide chemi- 
cal feedback between the fir st and second generations, and 
ISchaerer & Charbonnell(l201 ll) explore the effects of varying the 
slope of the high-mass end of the stellar initial mass function 
(IMF) on the resulting second-generation population. Requiring 
that the present-day ratio of first- to second-generation stars be 
1 :2, the authors calculate that initial globular cluster masses must 
have been 8 to 10 times their present-day values for a Salpeter 
IMF. They then compare the mass in first-generation stars that 
was lost from globular clusters at early times to the mass of 
the present-day Galactic stellar halo. Since roughly 35% of the 
mass in a Salpeter IMF is in low-mass, long-lived stars, and the 
present-day Galactic globular cluster system comprises 2% of 
the mass of the halo, they estimate that escaped first-generation 
globu lar cluster stars make up 5 - 8% of the mass of the halo. 

ISchaerer & Charbonnell d201 ll) then consider the origin of 
the second-generatio n stars found in the halo by MG10 and by 
ICarretta et al.l d2010l) . specifically whether they can reasonably 
be understood as having been lost to the halo during the strong 
early mass loss that removed roughly 90% of the first-generation 
stars or whether they were transferred during later episodes of 
globular cluster tidal disruption. Folding in some assumptions 
involving the initial cluster mass function (ICMF), they con- 
clude that all second-generation stars in the halo were lost at 
early times. This raises the fraction of cluster mass that was lost 
in this early phase, which necessarily increases the initial cluster 
mass to 15 to 25 times the present-day mass and also raises the 
fraction of halo stars that originally formed in globular clusters 
to as much as 20%. 

However, given the observed examples of globular cluster 
dissolution in progress, e.g., Palomar 5 (Rockosi et al. 2002 ; 
Odenkirchen et al. 2003), N GC 5466 (|Belokurov et"aT1 l2006h 
and the GD-1 stellar stream dGrillmair & Dionatosl 120061) . as 
wel l as the relative emptiness of the more inhospitable regions of 
the Gnedin & Ostrikerl d 19971) ' 'vital diagram", we expect clus- 
ter dissolution to be the dominant source of halo field stars with 
second-generation abundance patterns. Theoretical studies of the 
ICMF, early star cluster mass loss, and the evolution of the clus- 



ter mass function (e.g., Baumgardt et al. 2008; Parmentier et al. 
2009; Decressin et al. 2010) all suggest that far more star clusters 
were initially formed than have survived to the pr esent day, espe- 
cially at the lower end of the ICMF. The study of Jordi & Grebell 
(2010), which used matched-filter techniques to search for faint 
tidal tails around Galactic globular clusters in the SDSS imaging 
footprint, found stars with cluster-like photometry outside the 
tidal radii of a number of clusters without pronounced tidal tails, 
indicating that stars continue to escape from globular clusters 
even in the absence of dramatic tidal features. In the limit where 
cluster dissolution is the sole source of second-generation field 
stars, we can estimate the number of present-day globular clus- 
ters that would need to be disrupted to provide the observed 2.5% 
of halo stars with second-generation chemistry, Nd. For a typi- 
cal present-day 1 : 1 ratio of first- to second-generation stars, and 
assuming a stellar halo mass of 1O 9 M (e.g., Freeman & Bland- 
Hawthorn 2002) and a typical present-day globular cluster mass 



of 5 x 1O 5 M , N d =* 



2xf. 2G xM hl ,i„ 



- 100. Since the current number 



l0 ' i,a_ <M gc > 

of of globular clusters is on the order of 150, this means that the 
initial population of globular cluster s was significantly larger, 
but th is is not an outrageous claim: iMackev & van den B ergh 
(120051) calculate that the ratio between the initial and present- 
day number of globular clusters in the Galaxy is 3/2, whereas 
we derive 5/3. 



Whatever the details, it seems clear that first-generation stars 
lost from globular clusters at early times are an important contri- 
bution to the construction of the Galactic halo, if two-generation 
self-enrichment scenarios for globular cluster formation are in - 
deed the correct model. The models of Decress in et alJ d2007l) 
D'Ercole et al. (2008), Carretta et a i1 d2010l) . and lVesperini et al.l 



Q2010) all predict that 90% of the first gener ation ought to b e lost 
from globular clusters at early times, while Conrov (201 1) esti- 
mates that as much as 95% of the first generation ought to be lost. 
The mass of first-generation stars lost during this phase by the 
Ngc globular clusters in the present-day system can be expressed 

1). Considering the early 



as Mj(? 



M GC system x 5 x ( 



i-fk, 



contributions of first-generation stars from globular clusters that 
have survived to the present day as well as the clusters that have 
completely dissolved, we estimate that f ° c has a minimum value 

of + x m£) + M cc sy Stem = 2% X (| + f ) = 17% 
of halo field stars, with both first- and second-generation abun- 
dance patterns, originally formed within globular clusters, plus 
an additional unknown contribution of first-generation stars from 
clusters too low-mass to self-enrich that have dispersed by the 
present day. 
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5. Future challenges 

Based on the current two-generation models for globular clus- 
ter formation, the early phase of globular cluster formation must 
have contributed a significant number of first-generation stars to 
the Galactic halo. From the data presented herein, by MG10, and 
bv ICarretta et al. (2010), it is also clear that second-generation 
globular cluster stars are a component of the halo field. Further 
investigations of the process of globular cluster formation will 
help to clarify the role that globular clusters played in the early 
assembly of the Galactic halo. There are two specific direc- 
tions that new observations could take that would be particularly 
insightful: investigating the relationship between the minimum 
mass for cluster self-enrichment and the larger galactic environ- 
ment, and searching for star clusters still in the process of self- 
enrichment. 

These studies will necessarily require observations of star 
clusters outside the Milky Way, as a way to compare cluster 
formation in different large-scale environments. Observing very 
young star clusters may offer the opportunity to see the process 
of self-enrichment in progress. Carretta et alJd2010t) suggest that 
clusters in the process of violent relaxation (at a few xlO 7 years 
old) should have a central compact cluster, an extended, non- 
bound halo and outflowing gas. It is suggested in Conroy (201 1) 
that th e extended clusters discovered in M3 1 by iHuxor et aH 
J200l are in this state (see also Huxor et al. 2011). If this is 
the case, the stars at large radius in the extended clusters should 
have almost entirely first-generation chemistry, and strong radial 
gradients in light-element abundances, most readily observable 
in CN and CH molecular features, should be present. 

To explore the question of the minimum cluster mass for 
self-enrichment, it w ould be useful to c o nsider star clusters in 
lower-mass galaxies. IConrov & SpergeJ d201 ll) claim that the 
weaker tidal field of the Large Magellanic Cloud (LMC) ought 
to permit the formation of a second stellar generation at lower 
cluster mass than in the Milky Way. Spectroscopic observa- 
tions of red giant-branch stars in intermediate-age populous clus- 
ters in the LMC would be particularly interesting because re- 
cent photometric work (e.g., Mackey et al. 2008; Milone et al. 
2009; Goudfrooij et al. 2011) has uncovered broadened or split 
main-sequence turnoffs in a large fraction of them. Since they 
are much younger than Galactic globular clusters (typically 1-4 
Gyr), age differences on the order of a few hundred million years 
are visible at the turnoff, and this age difference is suggestively 
similar to the age difference expected between first- and second- 
generation stars in Galactic globular clusters. Obtaining spectra 
of turnoff stars at the distance of the LMC would be extremely 
difficult - indeed, high-resolution spectroscopic studies of RGB 
stars in these clusters are limited to fairly small samples (e.g., 
Mucciarelli et al. 2008). However, a lower-resolution study using 
a multiobject spectrograph could collect CN and CH data analo- 
gous to the data presented in this paper, and would allow a search 
for star-to-star variations in carbon and nitrogen abundance for 
a larger data set per cluster. A modified version of the SSPP has 
been developed for application to spectra with resolving power 
as low as R ~ 1000), and is already being used with a variety of 
low-resolution data from non-SDSS sources (e.g., see Li et al. 
2010; Humphreys et al 2011). These data would also permit an 
investigation of theoretical claims (e.g., Conroy & Spergel 201 1) 
that lower-mass galaxies permit the formation of multiple stellar 
generations in lower-mass star clusters, by comparing the CN- 
CH behavior of intermediate-age LMC clusters across a range 
of masses. 
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